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Abstract. In spin crossover materials, an abrupt phase transition between a low 
spin state and a high spin state can be driven by temperature, pressure or by light 
irradiation. Of a special relevance are Fe(II) based coordination polymers where, in 
contrast to molecular systems, the phase transition between a spin S = and a spin 
S = 2 state shows a pronounced hysteresis which is desirable for technical appHcations. 
A satisfactory microscopic explanation of this large cooperative phenomenon has been 
sought for a long time. The lack of X-ray data has been one of the reasons for the 
absence of microscopic studies. In this work, we present an efficient route to prepare 
reliable model structures and within an ab initio density functional theory analysis 
and effective model considerations we show that in polymeric spin crossover compounds 
magnetic exchange between high spin Fe(II) centers is as important as elastic couplings 
for understanding the phase transition. We discuss the relevance of these interactions 
for the cooperative behavior in these materials. 
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1. Introduction 

An intensively debated class of materials with potential applications as optical switches, 
sensors or memory devices [HEllSlll] are spin-crossover polymer (SCP) systems involving 
transition metal ions linked with organic ligands [5]. These systems show a sharp 
transition triggered by variation of temperature, pressure or by light irradiation between 
a low-spin (LS) ground state and a high-spin (HS) excited state with a thermal hysteresis 
loop [6j. Especially important in these materials is the large cooperativity shown at the 
HS-LS transition in contrast to molecular spin crossover systems, for example. The 
origin of this transition and its cooperativity has been mainly discussed in the frame of 
elastic models [6l [7l [H [H [10] , and only recently a possible role of magnetic exchange was 
suggested [TT], [T2]. Still, a conclusive ab initio microscopic study where all important 
interactions are considered is missing and the origin of the large cooperativity has not 
been completely settled. It is our purpose to investigate this issue in what follows. 
Ab initio theoretical studies for SCP systems are faced with major difficulties due to 
the nonexistence of well-determined crystal structures. To our knowledge, electronic 
structure calculations have only been performed at the level of semiempirical extended 
Hiickel approximation for an idealized triazole-bridged Fe(II) chain [I3]. In the present 
work, we overcome the unavailability of structural data by predicting a crystal structure 
for a Fe(II) spin-crossover polymeric crystal using known experimental constraints and 
a combination of classical force field and quantum mechanical molecular dynamical 
methods. We analyze with Density Functional Theory (DFT) calculations the LS-HS 
phase transition and show that there exists an interplay between magnetic exchange 
and elastic properties that is responsible for the large cooperativity in these systems. 
We also corroborate the quality of our designed structure by comparing with magnetic 
experiments done on a real sample. Our methodology and results provide a new 
perspective on the parameters underlying the traditional theoretical approaches. 

There have been a number of attempts to theoretically account for the features 
of the HS-LS transition in spin crossover materials. Most of the theoretical work is 
based on elastic considerations. Two types of elastic models that focus on the faithful 
reproduction of macroscopic quantities like the HS fraction, are the following: In the 
first approach, the cooperativity in the HS-LS transition is defined in terms of local 
distortions which interact with one another elastically causing a long range effective 
interaction between spin states. This leads to an Ising-type Hamiltonian [HI [15], \16\ 
H = J2i j JijO'iCTj which describes the elastic interaction between spin states (LS and 
HS) in terms of fictitious spin operators (cr = -1 (1) for LS (HS)) coupled via a nearest 
neighbor interaction Jij. The coupling constants Jij are parameters of the theory and 
have not yet been determined from a microscopic model. In an alternative approach, the 
free energy of spin crossover systems is calculated based on an anisotropic sphere model 
that describes volume and shape changes of the lattice at the transition [7l[8l[9]. One of 
the few attempts to include magnetic interactions is the recent proposal done in Ref. [I2] 
where the author considers a one-dimensional model for HS-LS systems which contains 
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elastic and magnetic Ising exchange interactions. The ground-state phase diagram is 
then obtained by the transfer matrix technique for different relative elastic to magnetic 
coupling strengths. 

In the present work we concentrate on the electronic and magnetic degrees of 
freedom in a spin crossover polymer, and investigate their influence on the microscopic 
origin of the cooperativity in the HS-LS transition. While it has been assumed in the 
past that the Fe(II) center nearest neighbor interaction is mostly of phononic origin, our 
study indicates that a signiflcant part of this interaction arises from magnetic exchange. 

2. Crystal Structure 

For this investigation it is indispensable to obtain a reliable crystal structure 
suitable for DFT analysis. We aim at describing the complex Fe[(hyetrz)3](4- 
chlorophenylsulfonate)2 ■ 3H2O [I7] {hyetrz stands for 4 — (2' — hydroxy ethyl) — 1,2,4 — 
triazole) (see compound 1 in Fig.[l]), which was synthesized from 2-hydroxyethyltriazole 
and iron(II)-p-chlorobenzenesulfonate as described in Ref. [T8] . 

This compound precipitates as a flue, polymeric powder and single crystals cannot 
be grown because the polymer is insoluble in water and organic solvents [H]. Melting 
or sublimation attempts result in decomposition. Hence, X-ray structure analysis is not 
possible and even the X-ray powder diagram consists of a few broad peaks only which 
prevent the structure to be determined from X-ray powder data. Other polymeric Fe- 
triazole compounds have similar properties and no single crystal structure for these 
systems is known. In the literature only single crystal structures for trimeric Fe 
compounds (e.g. Ref. [201), and polymeric Cu-triazoles [211 [22] were found. In these 
compounds, the metal ion is coordinated by six nitrogen atoms and neighboring metal 
ions are connected by three pyrazole bridges. For the polymeric Fe triazoles a similar 
structure is assumed, see Fig. [H This structure is also supported by spectroscopic 
methods [22] including solid state NMR [23] . 

In view of the above and the requirement of having reliable crystal structures 
for microscopic studies, we design on the computer a model system of polymeric Fe 
triazole, as a basis for calculating the electronic and magnetic properties. We consider all 
experimental information available and construct a crystal structure as close as possible 
to the actual structure. We employ a method that has been developed and tested on 
the coordination polymer Cu(II)-2,5-bis(pyrazole-l-yl)-l,4-dihydroxybenzene which has 
a simpler structure and less atoms per unit cell than polymeric Fe-triazole [24l [25]. 

Since we aim at understanding the HS-LS transition with accurate all-electron DFT 
calculations, which are computer intensive, we keep the essential features of the material 
and simplify those elements that are secondary to the transition, like the nature of the 
substituent R and the counter ion X~ (spin transitions are observed for a wide range of 
different substituents R and counter ions X~). We consider compound 2 in Fig. [U with R 
= CH3 and X~ = F~. This model structure has the short range environment of the Fe(II) 
centers exactly as in structure 1 in Fig. [U while the longer range environment (> 5 A) 
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Figure 1. Chemical diagram of polymeric Fe(II)-triazole. The model 
compound 2 differs from the synthesized compound 1 in the simplified R and X~ 
groups. 



of the Fe^"*" centers is significantly simplified. The molecular geometry of compound 2 
was constructed according to data from single crystal analysis of dimeric and trimeric 
Fe-triazole compounds p6]. A hypothetical, but sensible crystal structure of compound 
2 was built up with the minimum number of atoms per unit cell (72 atoms). The iron- 
triazole chain itself has Q^/m2/m2/ c symmetry. A crystal structure with hexagonal 
symmetry is in principle possible, e.g. in space group PG^/m, but we chose a structure 
with P2i/m symmetry with two formula units per unit cell. A similar arrangement of 
chains is also found in a corresponding Cu polymer, [Cu(hyetrz)3](CF3S03)2-H20 [21] . 
(This structure is triclinic, with space group PI, Z = 2, but the deviations from 
monoclinic symmetry are probably caused only by the symmetry-breaking hydroxyethyl 
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ligands and CF3SO3 counterions. Otherwise the symmetry would be P2i/m.) In our 
structure, the Fe^"*" ions are located on crystallographic inversion centers, whereas all 
triazole units contain a crystallographic mirror plane between the neighboring nitrogen 
atoms. All Fe^"*" ions are crystallographically equivalent and we enforce, for simplicity, 
a perfectly octahedral environment of the Fe^^ ions. 

In the design and analysis of the model structures we employ four distinct classical 
and ab initio methods. For preoptimization of the model structures we consider 
a classical modified Dreiding force field [27] with atomic charges calculated by the 
Gasteiger method [2H| as implemented in the Cerius^ package [29j. We then perform 
quantum mechanical first principles calculations within DFT with three different basis 
sets, each for a distinct purpose: The plane wave ab initio molecular dynamics (AIMD) 
method [30] is used for determination of precise equilibrium structures. The linearized 
augmented plane wave (LAPW) method [31] is used to determine accurate electronic and 
magnetic properties and the linearized / Nth order muffin tin orbital (LMTO/NMTO) 
methods [32] are used to calculate effective Fe Wannier functions and to understand the 
low energy excitations of the system. 

For the AIMD calculation [30] we considered a plane wave cutoff of 30 Ryd for the 
plane wave part of the wave function and we used the following sets of (s,p,d) projector 
functions per angular momentum: Fe(2,2,2), F(2,2,l), N(2,2,l), C(2,2,l) and H(2,0,0). 
We employed a (4 x 4 x 4) mesh and the P21/m symmetry was preserved during the 
relaxation with the help of 131 constraints. 

For the LAPW calculations [3l], we employed a (6 x 6 x 9) mesh and a value 
of i?MT-f^max = 4.2 that is sufficiently large due to the small radii Rmt of the hydrogen 
atoms. 

The NMTO-downfolding [32] calculations, which rely on the self-consistent 
potentials derived out of LMTO [33] calculations, were carried out with 40 different 
empty spheres in addition to atomic spheres to space fill. The convergence of LMTO 
calculations in each case was cross-checked with full potential LAPW calculations. 

All DFT calculations were performed within the Generalized Gradient Approxima- 
tion [34]. 

We first proceed with the relaxation of the model structure with the classical force 
field within the P 1 1 2i/m space group (non-standard setting of P2i/m), keeping fixed 
the Fe-N distances of the FeNe octahedra to the values (ipe-N ^ {2.00 A, 2.05 A, 2.08 A, 
2.09 A, 2.10 A, 2.12 A, 2.15 A, 2.20 A}. This constraint determines the value of crystal 
field splitting from the onset and constitutes our control parameter as shown below. As 
the N-N bond of the triazole (see Fig. [1]) has experimentally a well defined bond length 
of (iN-N = 1-38 A, the choice of a Fe-N distance leads automatically to lattice parameters 
linear in (ipe-N- The relative change of the lattice parameter c (the chain direction) is 
larger than that for a and b. This means that the volume of the unit cell also increases 
linearly from V = 6508 A^ at dpe-N = 2.00 k to V = 7317 A^ at dpe-^ = 2.20 A. It 
should be noted that the force field we use is not suitable for relaxing the F~ anion into an 
equilibrium position due to the lack of appropriate parameters for F~ anions. Therefore, 
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the F~ counter-ion is put into a likely position and the optimization is left to the next 
step where we perform a precise relaxation of the structure with the help of ab initio 
molecular dynamics. This AIMD step is essential as we find that the force field relaxed 
structures still show no or very bad convergence in LAPW, indicating inappropriate 
positions of at least some atoms. We relax all unconstrained atoms within AIMD until 
the forces that initially are of the order of 1000 mRyd/Bohr until they are far below 
1 mRyd/Bohr. More significantly, we converge bond lengths to within 10~^ A and bond 
angles to within 0.1°. We cross check the final AIMD relaxed structures by calculating 
the LAPW forces and making sure that they are of the order of 10 mRyd/Bohr or less 
(with the only exception of the difficult to place counter ions F~ for which we allow 
30 mRyd/Bohr as they do not play a crucial role for the Fe(II) center interactions). 
The stipulation that forces are very low in two fundamentally different DFT codes 
is a high standard and gives us confidence that our conclusions about the electronic 
and magnetic structure are drawn on a solid basis. The precision to which we converge 
bonds and angles means that we can then proceed to predict interactions between Fe(II) 
centers with the necessary precision of a few Kelvin. In Fig. [2] a representative of the 
resulting structures is shown. The top panel pictures the chain of FeNg octahedra with 
alternating orientations, and in the bottom panel we demonstrate the arrangement of 
the Fe(II) chains in the crystal. 

3. Energy scales 

At the atomic level, two sets of energy scales are responsible for the LS and HS state 
of the Fe(II) centers, namely the crystal-field splitting and the Hund's exchange and 
Coulomb interactions. 

The crystal-field splitting, as mentioned above, is being fixed beforehand in the 
construction of the structures with given (ipe-N distances. In Fig. [3] we show the 
dependence of the crystal field splitting on the size of the FeNg octahedra calculated 
with LAPW (magenta symbols) and NMTO (green symbols). Note the good agreement 
between the two calculations. The crystal field splitting values are obtained by 
determining the first moment of both t2g and Cg densities of states from non-spin 
polarized LAPW calculations and through construction of an Fe d only Hamiltonian 
in case of NMTO calculations. The atomic state diagrams for Fe 3d schematically 
demonstrate the relationship between crystal field splitting A and the spin state {S = 
ov S = 2). The Hund's exchange is taken into account to some extent within the spin 
polarized-GGA approach. 

In the extended system, two type of interactions contribute to the phase transition; 
the phononic excitations and the exchange interaction J, due to nearest neighbor 
superexchange between Fe(II) centers which is typically antiferromagnetic. 

The competition between all these energy scales determines the nature of the phase 
transition and its cooperativity, as we will discuss in section [6l The role played by the 
phonons in driving the LS-HS transition in the spin-crossover systems has been discussed 
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Figure 2. Simplified Fe(II) triazole structure. Top: Individual chain of 
octaiiedrally coordinated Fe^+ ions. Note tlie strong linkage between Fe(II) centers 
via triple N-N bridges. The orientation of the FeNg octahedra is alternating along the 
chain. Bottom: Arrangement of polymeric Fe(II) chains in the crystal. 
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Figure 3. Crystal field splitting A as function of Fe-N distances. The results 
obtained with LAPW (magenta) and LMTO (green) agree very well. The pink (white) 
background symbolizes the respective colors of the LS (HS) compounds and indicates 
the A value region that leads to LS (HS). The atomic Fe 3d state diagrams visuaHze 
the relationship between crystal field splitting and spin state. 



at length in terms of elastic models [Tj [HI [IHl [HI [15]. In our model calculations, the 
phononic degrees of freedom are frozen and we investigate the role of the electronic and 
magnetic degrees of freedom. 

4. Electronic structure 



We performed spin resolved electronic structure calculations for these systems within 
DFT by considering the GGA [3l] as exchange-correlation functional and the LAPW 
basis. In Fig. [4] we present for four of the designed structures the projection of the 
total density of states (DOS) on the Fe 3d orbitals, which are responsible for the 
magnetism of the material. The first and second panels correspond to the structures 
with dpe-N = 2.00 A and rfpe-N = 2.08 A respectively and show a perfect symmetry 
between spin up (red) and spin down (blue) density of states and therefore define a 
LS state (5* = 0). The occupied states can be identified as the six t2g states, while the 
empty states are the four Cg states. The splitting A between t2g and Cg states diminishes 
with increasing (ipe-N from A = 2.66 eV to A = 2.22 eV, respectively (see also Fig. [3|). 
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Figure 4. Spin resolved density of states for selected Fe(II) triazole 
structures. Red and blue colors denote spin up and spin down species respectively. 
The two upper panels for low Fe-N distances show symmetric density of states 
corresponding to a low spin {S — 0) state. The two lower panels for high Fe-N distances 
with strong spin up spin down asymmetry correspond to a high spin {S — 2) state. 



The DOS behavior completely changes as the Fe-N distance increases to (ipe-N = 2.10 A 
(see the third panel of Fig. HI). Now spin up (red) t2g and Cg states are completely 
filled, and spin down (blue) t2g states show only a partial filling with one electron. The 
imbalance between up and down electron numbers is — = A which corresponds to 
the HS state (S* = 2). This situation remains if the Fe-N distance increases further to 
(ipe-N = 2.20 A, only the splitting A between t2g and Cg diminishes. 

Thus, by carefully preparing a series of model structures that correspond to the LS 
and HS sides of the spin crossover transition we manage to microscopically describe the 
LS-HS transition which occurs between the structures with rfpe-N = 2.08 A (S* = 0) and 
cipe-N = 2.10 A (5* = 2). 

In order to quantify energetically the HS-LS spin transition, we show in Fig. [5] 
the total LAPW electronic energies obtained within the spin-polarized GGA (sp-GGA) 
approach. We note that there is a discontinuous jump between the LS (S=0) energies 
and the HS (S=2) energies. The relative electronic energy differences between the 

HS and LS systems — — jTfgf — is about 10~ which agrees with the relative energy 
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Figure 5. Ground state energies for the set of model structures obtained 
within sp-GGA with the LAPW basis . 

estimates for spin crossover molecular systems [35]. Since the LS-HS phase transition 
occurs between structures (ipe-N = 2.08 A and rfpe-N = 2.10 A, we designed one more 
structure with (ipe-N = 2.09 A in order to probe the sharpness of the transition. While 
there are indications that this structure might represent an intermediate magnetic state 
with S ~ 1.5 per Fe(II) center, we do not include it in Fig. [5] as it is very hard to 
converge. This result indicates that the spin crossover transition in polymer systems 
occurs in a very narrow range of crystal field splittings, i.e. is very sharp, as observed 
experimentally [5j. 

5. Exchange interaction and magnetic properties 

For the magnetic behavior of this system, we derive first from the DFT electronic 
calculations a Hamiltonian which describes the effective interaction between Fe(II) 
centers. The NMTO-based downfolding method, designed to pick up selectively the 
low-energy bands from the whole group of LDA/GGA bands of a compound, has been 
used to construct the Fe 3d only Hamiltonian of the Fe-triazole compounds. This 
method has proven to provide reliable information on the low energy properties of 
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inorganic [36l [37l [38] and organic [25] transition metal compounds. The tight-binding 
basis in which these Hamiltonians are constructed form the set of "effective" functions, 
which span the Hilbert space of the Wannier functions corresponding to the low-energy 
bands. Fig. [6] shows the plot of one member of such a set, namely the downfolded 
Wannier function corresponding to Fe 3dxy In the figure, two such Fe 3dxy Wannier 
functions have been placed at two neighboring Fe sites. While the central part of such 
an effective function has the Fe Sd^y symmetry, the tails of the function are shaped 
according to integrated out degrees of freedom in the system, like C sp, N sp, F sp 
and H s. As is evident from the plot, substantial weight of these tails resides on 
neighboring triazole rings. The presence of these tails hints to an enhanced electronic 
hybridization between the adjacent Fe^+ ions, contributing to the cooperative nature of 
the HS-LS transition. Out of these calculations we can estimate the various hopping 
matrix elements, t, between the d orbitals of adjacent Fe(II) centers. The values of these 
hopping parameters range between 1 meV to 80 meV quantifying the strength of the 
various interaction paths between neighboring d Fe(II) orbitals. 

In order to get estimates of the magnetic superexchange coupling constants J 
between neighboring Fe(II) centers for the HS Fe-triazole structures, we considered two 
approaches: 

i) Total energies. We calculated within sp-GGA total energies of ferromagnetic 
and antiferromagnetic Fe^"*" spin configurations. Considering a spin-Hamiltonian 
H = JSiSi+i between nearest neighbors Fe^+ spins Si, the ferromagnetic and 
antiferromagnetic energies for two Fe^"*" ions in the unit cell of the Fe-triazole are 
given by E^m = 8J and Eafm = — 8J. Within sp-GGA Epu — -E'afm = 33 meV 
for the (ipe-N = 2.20 A structure and therefore J f» 2.1 meV = 24 K. ii) Perturbation 
theory |39] where J can be obtained in terms of the hopping parameters t between Fe(II) 
centers and the onsite Coulomb repulsion f/ as J ~ At'^/U. For the HS (ipe-N = 2.20 A 
structure, the significant t obtained within the downfolding method is 48 meV, and for 
U = 4-5 eV this gives J ~ 22K which is very similar to the value obtained with the 
difference of total energies. 

The results of our model calculations can be now compared with the magnetic prop- 
erties obtained on the real samples of polymeric Fe[(hyetrz)3](4-chlorophenylsulfonate)2- 
3H2O. Variable-temperature magnetic susceptibility measurements in the temperature 
range 2 — 350 K and magnetic fields 0.02 — 0.2 T were carried out on powder samples 
of Fe(II) triazole using a Quantum Design SQUID magnetometer MPMS-XL. In Fig. [7] 
we show the susceptibility x measurements where Tx has been plotted versus T. Our 
sample shows hysteresis at T = 80 K with a width of 20 K. Since this system consists 
of spin S = 2 Fe(II) chains with weak interchain interactions, we have analyzed the 
magnetic susceptibility in the frame of a spin S = 2 Heisenberg chain model. 

Note that the Fe(II) triazole SCP systems are Haldane S = 2 chains. A Haldane 
gap is expected to exist between the ground state and first excited states. The reason for 
not observing this gap, is that the transition temperature at which the HS-LS transition 
happens is higher in energy than the Haldane gap energy and therefore, the system goes 




Figure 6. Wannier functions for the Fe-Fe interactions. The chosen example 
shows two Wannier functions with 3dxy symmetry on two neighboring Fe(II) centers. 
The effective interaction (superexchange coupHng) between Fe{II) centers depends on 
the degree and relative sign of the overlap of the Wannier function tails on the pyrazole 
ring. 



into the non-magnetic 5 = phase before the gap in the S 
The susceptibility of an L-site chain is given by: 



XL 



2 chain can be observed. 



(1) 



k^T Tr [e-/^^] 

where H is the Heisenberg Hamiltonian H = JS'jS'j+i, fiB is the Bohr magneton, g is 
the gyromagnetic factor, /cb is the Boltzmann constant, T is temperature and S'f is the 
^-component of the spin on site i. In the thermodynamic limit, the bulk susceptibility 
at high temperatures can be obtained as a series expansion in 1/T: 
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For k-QT / J > 5(5* + 1) Eq. [2] compares very well to QMC data for spin 5 = 2 chains [40]. 
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Figure 7. Measured susceptibility. The solid line corresponds to a high-T fit 
as explained in the text. The susceptibility measured while lowering (increasing) the 
temperature shows a hysteresis of 20 K at the transition temperature of 80 K. 

The fit of Eq. (I2|) to the measured susceptibility of Fig. [7] is best for g = 2.2 and 
J = 11 K. 

The J values obtained from our ab initio calculations (J ~ 24:K) are slightly larger 
than the J extracted from the susceptibility data (J = 11 K) on the real sample, 
but remain in the same order of magnitude. Considering that i) we performed the 
calculations in a model structure and ii) the experimental measurements are affected by 
the quality of the samples, we can conclude that the comparison is quite good and the 
designed structures are reliable. 

6. Discussion 

One important topic of this work is the analysis of the various energy scales that 
contribute to the cooperativity of the HS-LS transition in SCP systems. In our 
calculations we froze the elastic degrees of freedom and concentrated on the electronic 
and magnetic properties for which we have quantified the corresponding parameters. 
A comparison with the elastic coupling constants estimated from Ising-like models [15] 
(in such an approach, a material with a transition temperature of T = 80 K would be 
described by elastic interactions of JeiasUc ~ 20 — 30 K) shows that in one-dimensional 
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Fe(II) triazole, the elastic coupling constants JeiasUc are of equal importance as the 
magnetic exchange for explaining the spin crossover transition, and the cooperativity 
should be understood as an interplay between elastic properties and magnetic exchange. 
When we cool the system from the HS state toward the HS-LS transition temperature, 
the elastic coupling tends to drive the system to the LS state, while the magnetic 
exchange tends to keep up the magnetic state for a larger temperature range (see Fig. [7]). 
In comparison, in the heating process, the magnetic exchange is initially absent (LS) 
and therefore the elastic interaction (vibrational phonons) initially drives the transition 
which has its on-set at a higher temperature than in the cooling process. The width 
of the transition between cooling and heating (hysteresis) is therefore enhanced by the 
magnetic interaction. 

A fundamental difference between the polymeric systems we are dealing with in 
this work, and molecular bi- (tri-, tetra-, . . . ) nuclear Fe systems is the connectivity 
between the Fe(II) centers. While the molecular systems [41] form isolated clusters 
of Fe(II) centers and therefore there is no strong connectivity between clusters, the 
polymers have important nearest neighbor interactions in the thermodynamic limit. 
This implies that for the polymers, the magnetic superexchange is not restricted to the 
cluster as in the molecular systems, but rather becomes important for the nature of the 
HS-LS phase transition. In Ref. [l2] various estimates of the magnetic J for molecular 
systems have been given. The values range between 4-6 K. The values we estimated for 
the present polymers are larger, between 11-24 K and in the energy range of the elastic 
constants, which indicates that the cooperativity in these SCP systems is most likely 
significantly enhanced by the exchange interactions. 

In conclusion, this work presents an efficient route to prepare reliable model 
structures for microscopic investigations and provides a new interpretation about the 
origin of the parameters underlying traditional theoretical approaches for the polymeric 
spin-crossover materials. The next step of this study will be the inclusion of lattice 
dynamics (phonons) in the ab initio calculations which is planned in our future work. 
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